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ABSTRACT 

We study the distribution of cosmic voids and void galaxies using Sloan Digital Sky 
Su rvey Data Release 7 (SD SS DR7). Using the VoidFinder algorithm as described 
by iHovle fc VogelevI (|2002D . we identify 1054 statistically significant voids in the 
northern galactic hemisphere with radii > 10h~^ Mpc. The filling factor of voids 
in the sample volume is 62%. The largest void is just over 30h~^ Mpc in effective 
radius. The median effective radius is 17h~^ Mpc. The voids are found to be sig- 
nificantly underdense, with density contrast S < —0.85 at the edges of the voids. 
The radial density profiles of these voids are similar to predictions of dynamically 
distinct underdensities in gravitational theory. We find 8,046 galaxies brighter than 
Mr — —20.09 within the voids, accounting for 7% of the galaxies. We compare the 
results of VoidFinder on SDSS DR7 to mock catalogs generated from a SPH halo 
model simulation as well as other A-CDM simulations and find similar void fractions 
and void sizes in the data and simu lations. This catalog is made publicly available at 
http://www.physics.drexel.edu/~pan/VoidCatalogi for download. 
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1 INTRODUCTION 

Redshift surveys of galaxies reveal a rich variety of large- 
scale structures in the Universe: clusters that span a few 
megaparsecs in radius, connected by filaments stretching up 
to many tens of megaparsecs, which in turn envelop vast 
underdense voids with radii of tens o f megaparsecs. Th ese 
large scale structures are described bv lBond et al.l (|l996t ) as 
a Cosmic Web of material that reflects the initial density 
fluctuations of the early Universe. While historically most 
attention has been paid by astronomers to the dense clus- 
ters and filaments, it is the voids that fill most of the volume 
in the Universe. These underdense regions strongly influence 
the growth of large scale structure. The statistics and dy- 
namics of cosmic voids and the properties of the few objects 
found within them provide critical tests of models of struc- 
ture formation. 

Observations of voids in the galaxy distribution have 
progressed as the depth, areal coverage, an d sampling d en- 
sity of galaxy redshift surveys has improved. iRoodI (| 19881 ') re- 
views the paradigm shift that occured beginning in the mid- 
1970s as the focus shifted from the study of galaxy surface 
distributions to three-dimensional spatial distributions pro- 
vided by redshift sur veys, and the i mpact of this revolution 
on studies of voids. iJoeveer et all (| 19781 ') identifled super- 
clusters and voids in the distribution of galaxies and Abell 



clusters. Pencil beam surveys of the C oma/Abell 1367 super- 
cluster |Gregory fc ThompsonI [l97l ') indicated large voids. 
iKirshner et al.l 1 198ll ) discovered a void in the Bootes region 
of the sky that is 50 h'^ Mpc in diameter, several times 
larger than any previous ly observed. T he Center for Astro- 
physics Redshift Survey (jHuchra et al.1,1983) and in partic- 
ular its extension to tu b = 15.5 ( de Lapparent et al.lll98^ : 
iGeller fc Huchralll989ll revealed that the large-scale struc- 
ture of galaxies is dominated by large voids and the sharp 
fllaments and wa lls that surround them. The Southern Sky 
Reds hift Survey (|da Costa et al. 19881: Maurogordato et all 
Il992h found similar results. The lGiovanelli fc Havned (|l985f ) 
survey detailed the supercluster and void structure of the 
Perseus- Pisces region. The de e per Las Campanas Re dshift 
Survey (|Kirshner et al.l Il99ll : IShectman et al.l Il996l ') con- 
flrmed the ubiquity of voids in the large-scale distribu- 
tion of galaxies. Comparison of optically-selected galaxy 
surveys with redshift surveys of infrared selected galaxies 
t Strauss et al . 1992; Fisher ct al. 1995; Saunders ct al. 200^; 
Ijones fc et al.ll20 04') indicated that the same voids are found 
regardless of galaxy selection. The comple ted Two Degree 
Field Galaxy Redshift Survey (2dFGRS; ICoUess fc etaT 



(2001")) and Sloan Digital Sky Survey (SDSS: lYork fc et al 
(2000); Abazaiian ct al. (2009)) now allow the most com- 
plete view to date of the detailed structure of voids. 

A variety of methods have been used to compile 
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catalogs of voids in both observations of galaxies or 
clusters and in simulations (using dark matter particles 
or mock galaxy catalogs). Deta iled discussion o f man y 
of these methods is given by IColberg fc etU] (|2008l ). 
who compare void finding techniques. For the purpose 
of finding voids in redshift survey observations, methods 
that are applicable to th e dist ri bution of gala-xies in - 



clude iKa uffmann & Fairall' (l99l'l; ' El-Ad fc Piraiil 
Aikio fc Maehoenen (199&); Hovle fc VqgelevI 



1997 ) 



20021 ) 



Nevrinckl (120081 ): lAragon-Calvo et al.1 l|2010l ).^amples of 



applications of such methods to galaxy redshift surveys 
include analyses of t he Southern Sky Redshift Survey 
l|Pellegrini et al.l [1989'), the first slice of the Center for 
Astrophysics Redshift Survey (|Slezak et all 1 19931 ). as 
well as the full extens ion of the CfA Redshift Survey 
l|Hovle fc VogelevI I2OO2I), the IRAS 1.2Jy a nd Optical 
Redshift Surveys (|E1-Ad et al.1 Il997l : lEl-Ad fc P iranJ997 
200c), the Las Campanas Redshi ft Survey (iMiiller et al] 
2000r). the IRAS PSCz Survey (IHovle fc VogelevI l2002l: 



Plionis fc BaailakosI l20oj). t he 2dFGRS jHovle fc VogelevI 
2004 ICeccarelh et al.l 12006: 'Tikhonov"200d^, and proli: 



nary data from the SDSS (Tikhonov 2007; Foster fc Nclso^ 
I2OO9I ). 

The importance of cosmic voids as dynamically-distinct 
elements of large-scale struc t ure i s clearly established by 
theory (iHoffman fc Shaham 1982; Hausman et al.l 



Fillmore fc Goldreich. 1984: Ickc 198 4: BertschingeJ 
Blumenthal et ahT 199 2: Shctli fc v an de Wevgaeru 



1982 



1985 



2004 



Patiri et al.l I2OO6I : IPurlanetto fc Piranl bOOfiP . Linear 



theory predicts that the interior of the voids should 
reach a fiat plateau and the boundaries of the voids 
should be quite sharp. Simulations of structure forma- 



tion (e.g.. iRegos fc Gelled (Il99ll ): iDubinski etHI J 19931) 



Golberg et al.l 



Ivan de Wevgaert fc van Kampen l|l99j ): 

(2005D) demonstrate that large voids are caused by super- 
Hubble outflows that are nearly spherically symmetric out 
to near the edges of the voids. Tidal effects of clusters only 
become important for objects near the walls around voids. 
Simulations of the Cold Dark Matter model for structure 
formation indicate that the interiors of voids should 
include dark matter filaments and many low mass halos 
jMathis fc Whitdl2002l : iBenson et aLll2003l : ICottlober erahl 
I2OO3I ). Identifying these structures is an important test of 
this model. 

The properties of large voids in the distribu tion of 
:alaxi es may provide st r ong t ests of cosmology. iRvdenI 
19951 ): iRvden fc Melot^ (|l996l ) discuss the use of void 
shapes in redshift space as a cosmological test. More re- 
cent work examines voids as a pro be of dark e nergy 



cent work exammes voids as a pr o be 01 dark e nergy 
JPark fc Led l2007l: iLee fc ParkI l2009l : iBiswas et al] I2OI0I : 
iLavaux fc Wandeltll201ciF Comparison of voids at low and 
high re dshift may prov ide a strong test of the ACDM 
model l|Viel et al.l I2OO8I ). The abundance of cosmic voids 
is a critical probe for non-gau ssianity in the initial con- 
ditions for structure forma tion (|Kamionkowski et al.|[2009l : 
IChongchitnan fc Silkl [20T0I ). Beyond tests of the ACDM 
model, the properties of voids galaxies may even constrain 
alternative theories gravity (Hui ct al. 200^). 

Mapping the voids is important both for studying large- 
scale cosmic structure and because they are a unique astro- 
physical laboratory for studying galaxy formation. Gravi- 
tational clustering within a void proceeds as if in a very 



low density universe, in which structure formation occurs 
early and there is little interaction between galaxies, both 
because of the lower density and the faster local Hubble ex- 
pansion. Goldberg fc Vogeley (2004) show that the interior 
of a spherical void with 10% of the mean density in a fiat 
^mattfir = 0.3 h = 0.7 uuiverse evolves dynamically like an 
^m atter = . 02, = 0.48, h = 0.84 uuiverse. 

iPeeblej (|200ll ) describes the "void phenomenon" : galax- 
ies of all types appear to respect the same voids, in contrast 
to the prediction of CDM that low density regions should 
contain many low mass objects. iTikhonov fc KlvpinI (|2009l ) 
find, using comparison of voids and void galaxies in the local 
volume with high-resolution simulations, that the emptiness 
of voids is a problem for ACDM. 

Voids are expected to harbor many low mass halos that 
are the ideal breeding grounds for faint galaxies; if the low 
mass halos predicted by CDM harbor luminous galaxies, 
then they should be optically visible. Optical observations 
have n ot revealed a large population of fainter galaxies in 
voids jThuan et al.lll987l : iLindner et al] Il996l : iKuhn et al.l 
ll997l : |Popescu et al.ll 19971 ). although the luminosity function 
in vo id s is shifted by about one full magnitude (.Hovle et al.l 
l2005h . [Tinker fc Conrovl (|2009l ) contends the ACDM void 
phenomenon is due to a lack of understanding of assem- 
bly bias as galaxies form, but their model predicts a 5- 
magnitude shift in maximum galaxy luminosity. If void halos 
contain gas, but too few stars to be visible, then then their 
gas might be detected. To date, blind HI surveys have not 
detected such a popul ation of HI rich but optically dark 
galajdes (|Havnesl |2008| ). Nearby Lyman-a clouds detected 
along lines of sight toward bright quasars show a strong 
preference for inhabiting the voids, but most of these clouds 
seem to be associated with galaxy structures (Pan et al., in 
preparation) . 

In contrast to a picture in which star formation in void 
halos is suppressed, our analyses of void galaxies in the SDSS 
DR2 and DR4 samples indicate that void galaxies are bluer 
and have higher specific star formation rates than galaxies 
in denser environments (Hoi as et al. l l2004l2005l : |Park et al.l 
l2007i) . For the small number of dwarf galaxies in the earlier 
samples, we note even stronger trends with environment; at 
fixed morphology and luminosity, the faintest void galaxies 
are bluer and have higher star formation rates. Focusing on 
the blue population in voids, we flnd in these preliminary 
SDSS analyses, and Ivon Benda-Beckmann fc Mulled (l2008h 
find in 2dFGRS, that this blue population is not only more 
numerous, but also bluer and with higher star formation 
than in denser regions. 

In lHovle et al] (|2005l ) we flnd a much fainter exponen- 
tial cutoff in the luminosity function in voids {AM* — 1.1 
mag) but no evidence for a change in the faint end slope 
between voids and "walls." However, the uncertainties at 
faint magnitudes are quite large. We could not flnd a sub- 
population of "wall" galaxies selected by color, surface 
brightness profile, or Ha equivalent width that matched 
both the faint end slo pe a and charact eristic magnititude 
M* of void galaxies. In lPark et al.l (|2007t ) we again flnd that 
M* monotonically shifts fainter at lower density, and that 
the faint end (measured only down to Mr = —18.5) slope 
varies signiflcantly wi th density. These r esults are co nsistent 
with earlier analyses ICrogin fc GeUeij (|l999l . |2000| ). These 
trends also persist into the "wall" regions closest to large 
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voids ijCeccarelli et al.ll2008l ). When we estimate the mass 
function of void galaxies in SPSS and compare to the lumi- 
nosity function /Goldb erg et alJ boOSl) we find a good match 
with the predictions of the Press-Schechter model, thus the 
void galaxies appear to be nearly unbiased with respect to 
the mass. 

While we see some clear trends, controversy per- 
sists in the literature as to whether or not galaxies 
in voids differ in their internal properties from simi- 
lar objects i n denser regi o ns. F o r example, iRoias et alJ 
ll200 4'.'2005'):'Blanton et al.' ('20051): ' Patiri et al.l (|2006l ). and 
Ivon Bcnda-Bcckmann & Miillcr (200E) reach varying con- 
clusions that clearly depend on how environment is defined 
and which observed properties are compared. There is a 
marked difference between properties of the least dense 30% 
of galaxies (in regions with density contrast S < —0.5) and 
objects in deep voids which form the lowest density 10% of 
galaxies (in regions of density contrast 5 < —0.8, which is 
the theoretical prediction for the interiors of voids that are 
now going non-linear). All of these results, and the possible 
controversy among them, highlight the importance of build- 
ing the largest possible, publicly released catalog of voids 
and void galaxies. 

Lastly, for the purpose of examining the influence of 
environment on galaxy formation and evolution, it is impor- 
tant to make a distinction between void galaxies and iso- 
lated galaxies. Void galaxies are galaxies that reside within 
large scale void structures in the Universe. While this has an 
overall effect on the local environments of these void galax- 
ies, it does not preclude galaxies from residing within small 
scale dense environments, or clo ud-in-void as described in 
ISheth fc van de Wevgaer'3 (|2004l ). Isolated galaxies are gen- 
erally found by nearest neighbor distance measures typ- 
ically on the scale of small (Mpc) nearby e nvironments 
IIKarachentsev et al ]|2010l : lKarachentseTOlll973h ; they do not 
necessarily reside in large scale voids. 

The purpose of this paper is to obtain the largest void 
catalog available to date, for the purpose of allowing preci- 
sion cosmological tests with voids and more accurate tests 
of galaxy formation theories. We utilize a galaxy based v oid 
finding algorithm, "VoidFinder" (|lfovle fc Vogelevll200j '). to 
identify voids in the final galax;y catalog from SDSS (DR7). 
This void finding technique is shown to accurately identify 
large-scale cosmic voids with properties similar to those pre- 
dicted by gravitational instability theory. Section 2 describes 
the VoidFinder algorithm. Section 3 describes the SDSS data 
used for this research. Section 4 presents results on the var- 
ious properties of the voids found. Sections 5 and 6 describe 
several methods used to test the robustness of the method. 



2 VOIDFINDER 

VoidFinder is a galaxy-based void finding algorithm that 
uses redshift data to find statistic ally significant cosmic 
voids. VoidFinder is base d on the IEI- Ad fc PiranI (|l997l ') 
method as implemented bv llfovle fc VogelevI ( 20021 ) and uses 
a nearest neighbor algorithm on a volume limited galaxy 
catalog for void finding. This approach is highly effective in 
identifying large voids of density contrast S < —0.9 and ra- 
dius R > 10/i~^Mpc. The method is robust when applied to 
different surveys that cover the same volume of space (we 



have applied VoidFinder to IRAS PSCz, C fA2-HSSRS2, 2dF 
SPSS , 6dF and compared overlaps; see iHovle fc VogelevI 
(200^, [ 20041 )'). Our tests on cosmological simulations demon- 
strated that this method works in identifying voids in the 
distributions of bot h simulated galaxies and dark matter 
(|Benson et al.ll2003l ). 

VoidFinder is applied to volume limited galaxy samples. 
The galaxies are initially classified as wall or field galaxies. 
A field galaxy is a galaxy that may live in a void region 
whereas wall galaxies lie in the cosmic filaments and clus- 
ters. The distance parameter d for determining whether a 
galaxy is a wall or field galaxy is based on the third near- 
est neighbor distance (ds) and the standard deviation of the 
distance (fdg): 

d = ds + 1.5(7ti3 

In our galaxy sample, this selection parameter is d > 6.Sh~^ 
Mpc for field galaxies. With this value of d and choice of 
Mum = —20.09, all field galaxies reside in underdense re- 
gions with density contrast 5p/p < —0.47. Voids are ex- 
pected to be significantly underdense, containing approxi- 
mately 10% of the cosmic mean density. Near the edges of 
the voids, the density is expected to rise very sharply, dras- 
tically going from 20% of the mean density to 100%. Using 
this criterion for the edges of voids, it is expected that the 
distance criterion for void galaxies will depend on the density 
at the edge of the void and the spatial correlation of galaxies 
in voids and the fact that we are sitting on a galaxy. If we 
assume that the density at the edge of a void is 20% of the 
mean, then the expected density p around a galaxy near the 
void edge can be calculated as 



p(r)/p=(0.2)(e(r) + l) 



(1) 



where ^(r) is the two point autocorrelation function of the 
galaxy sample. The average density in a sphere of radius R 
around a galaxy near the void edge is therefore 

p(i?)/p= (0.2) (e(>) + l) (2) 

where ^{R) is the average value in a sphere of radius R. 
Over the scales of interest here, the redshift-space correla- 
tion function for galaxies in our volume-limited sample can 
be approximated by a power law, 



C = (s/soY 



(3) 



with So = 7.62 ± 0.67 and 7 = 1.69 ± 0.1. Using these values, 
we can determine the values of R20, the radial distance from 
a void galaxy where we would expect to encounter 20% of 
the mean density, and 5d=6.3, the expected underdensity of 
a void if its third nearest neighbor is found at a distance of 
6.3 Mpc away. We find 



R20 



+"-'^''h-'Mpc 



-0.74 



(4) 

-0.88 (5) 

We expect that at the edges of the voids 5 — —0.8, and in the 
centers of the voids 5 = —0.9. Thus, our choice for the value 
of d allows us to pick out void galaxies conservatively, se- 
lecting mostly galaxies that only live near the centers of the 
voids and not allowing void regions to grow into the nonlin- 
ear regime. All galaxies with third nearest neighbor distance 
da > 6.3/1-^ Mpc are considered to be potential void galax- 
ies and are removed from the galaxy sample, leaving us a 
list of wall galaxies. 
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We map out the void structure by finding empty spheres 
in the wall galaxy sample that remains. Wall galaxies are 
gridded up in cells of size 5 Mpc, which allows us to find 
all voids larger than 8.5 Mpc in radius. All empty cells 
are considered to be the centers of potential voids. A maxi- 
mal sphere is grown from each empty cell, but the center of 
the maximal sphere is not confined to the initial cell. Even- 
tually the sphere will be bound by 4 wall galaxies. There is 
redundancy in the finding of maximal spheres, but this is 
useful to define non-spherical voids. 

The sample of empty spheres now represents our poten- 
tial void regions. We sort the empty spheres by size starting 
with the largest. The largest empty sphere is the basis of the 
first void region. If there is an overlap of > 10% between an 
empty sphere and an already defined void then the empty 
sphere is considered to be a subregion of the void, otherwise 
the sphere becomes the basis of a new void. There is a cutofT 
of 10 Mpc for the minimum radius of a void region as we 
seek to find large scale structure voids that are dynamically 
distinct and not small pockets of empty space created by a 
sparse sample of galaxies. Any field galaxies that lie within a 
void region are now considered void galaxies. For further de- 
tails of the VoidFinder algorithm see iHovle fc VogelevI (|2002l . 
|2004| ). 



3 DATA: SDSS DR7 



We u se the SDSS Data Release 7 (DR7) (jAbazaiian et al] 
120091 ) sample of galaxies. The SDSS is a photometric and 
spectroscopic survey that covers 8,032 square degrees of the 
northern sky. Observations were carried out using the 2.5m 
telescope at Apache Point Observato ry in New Mexico in 
five photometric ban ds: u, g, r, i, and z iFukugita et al.ll 19961 : 
iGunn fc et al.|[T998l ). Follow up spectroscopy was carried out 
for galaxies with Petrosian r band magnitude r < 17.77 
after each photometric image was reduc ed, calibrated an d 
classified (|Lupton et al.ll200ll . [l99a : .Strauss fc et al.|[200^ '). 

Spectra were taken using circular fiber plugs with an 
angular size of approximately 55 arc seconds. If two galaxies 
were closer than this, we could only obtain the spectra of 
one; the other object is omitted unless there is plate overlap. 
iBlanton et al.l l|2003h addresses the issue of fiber collisions by 
assessing the relation between physical location of the galaxy 
and photometric and spectroscopic properties and assigns a 
redshift to the object missed by SDSS. 

We use the Korea Institute for Adva nced Study Value- 
Added Galaxy Catalog (KIAS-VAGC) (|Choi et aUlioiol) . 
Its main source is the New York University Value-Added 
Galaxy Catalo g (NYU-VAGC) La rge Scale Structure Sam- 
ple (brvoidO) l|Blanton et all [20051 which includes 583,946 
galaxies with 10 < r < 17.6. After removing 929 objects that 
were errors, mostly deblended outlying parts of large galax- 
ies, including 10,497 galaxies excluded by SDSS but that 
were part of UZC, PSCz, RC3, or 2dF, and also including 
114,303 galaxies with 17.6 < rur < Vim from NYU-VAGC 
(fuUO), there is a total of 707,817 galaxies. This catalog ofi'ers 
an extended magnitude range with high completeness from 
10 < r- < 17.6. There are 120,606 galaxies with z < 0.107 
and Mr < —20.09 in the volume limited sample. 




Figure 1. 10 /i ^ Mpc thick slab through the middle of the 
largest void at RA = 226.52960, DEC = 60.41244. Locations of 
galaxies (M^ < —20.09) are shown with *, and the locations of 
void galaxies are shown with -|-. The circles show the intersection 
of the maximal sphere of each void with the midplane of the slab. 



4 MEASUREMENT OF VOID PROPERTIES 

We identify 1,054 voids in SDSS DR7. The largest voids are 
30 Mpc in effective radius, where the volume of the void 
region is equal to the volume of the sphere with radius re//, 
and the median effective void radius is 17 Mpc. The 
voids cover 62% of the volume in the sample, and contain 
7% of the volume limited galaxies. We also identify 79,947 
void galaxies with SDSS spectra that lie within the voids in 
the r < 17.6 magnitude limited catalog, which corresponds 
to 11.3% of the magnitude limited galaxies. 



4.1 Void Sizes 

Figure [T] shows a redshift slab of 10 Mpc in thickness 
going through the center of the largest maximal sphere de- 
tected by VoidFinder. The intersections of the plane with all 
maximal spheres of void regions are shown. It can be seen 
that even with just the maximal spheres, a large volume of 
space is underdense and galaxies cluster strongly in large 
filament-like structures. Figure [2] (top) shows the radius his- 
togram based on the largest maximal sphere that defines the 
void region. It can be seen that the majority of the voids 
are small in size with a few very large void regions. Figured 
(bottom) shows the effective radius of the individual void re- 
gions. It is important to remember that the maximal spheres 
are limited to r > 10h~^ Mpc and thus only spherical void 
regions are found around 10 Mpc in effective radius. 
Most voids are not spherical and the skew in the effective 
radius histogram refiects the ellipticity of the voids. 

In figure |3] we see that the majority of volume occu- 
pied by voids are occupied by moderately sized voids. Even 
though a large number of voids are smaller in size, the actual 
volume distribution indicates that there is a preferred size 
for large scale structure in the Universe. This is consistent 
with observations sta rting with the CfA red shift survey to 
SDSS. As indicated in lShandarin et all l|2006l ). void sizes are 
largely determined by the cosmology. 



4.2 Radial Density Profiles 

The radial density profiles of the cosmic voids show that 
voids are significantly underdense, having less than 10% of 
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Figure 2. Distribution of void sizes as measured by the radius of 
the maximal enclosed sphere (top panel) and by effective radius 
(bottom panel). There is a cutoff of 10 Mpc for the holes that 
make up the voids and voids with rmax near this cutoff make up 
the majority of the void sample by number. The shift in the void 
distribution from the top to bottom panels indicates that the the 
void volumes are not well described by their maximal spheres; 
most voids are elliptical. Thus, the lack of small voids in the right 
panel is attributed to their ellipticity. 
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Figure 4. The average radial density profile of all 1,054 voids in 
the void catalog after scaling the profiles by Rvoid ^^<i stacking. 
The figure on the top is the profile of the enclosed volume, and 
the figure on the bottom is the profile in spherical shells. In both 
figures, there is a very sharp spike near the edges of the voids. 
The steep rise in the density contrast is because walls of voids 
are well defined. The peak at the edge of the void in the spherical 
shelss may be a feature of the density of the sample. 



maximai r {h'^ IVlpc) 



20 22 24 



errecTive r (Mpc n-i) 

Figure 3. Distribution of void sizes as a percentage of the volume 
occupied by the voids. As in Figure 2, the top panel sorts voids 
by their maximal sphere radii, on the bottom by their effective 
radii. Large voids occupy most of the volume with 50% of the 
volume occupied by voids with maximal sphere r > 13.8h~^ Mpc, 
and void size effective r > n.8h~^ Mpc. Note the peak of the 
radius histogram distribution around 22h~^ Mpc, the typical size 
of voids in the Universe. 



the average density all the way out to the very edge of the 
voids. Figure |4] (top) shows the stacked radial density pro- 
file of voids. The density is calculated from the volume en- 
closed to the given effective radius of the void. Figure |4] 
(bottom) shows a similar stacked radial density profile of 
the voids. However, the density is now calculated for spher- 
ical annuli. It can be seen that the walls of the voids are 
quite sharp, quickly growing from 10% of the average den- 
sity to 100%, and the voids are very well defined in terms of 
their density contrast with the outside Universe. It is clear 
then that these voids are distinct features of the Universe. 
A comparison with linear gravitatio n theory (Figure [51 re- 
produced from lSheth fc van de Weygaertl (| 20041 ')') shows the 
sam e "bucket shaped" rad i al den sity profile (see also Figure 
4 of iFiUmore fc GoldreichI (|l984l ')'). 



4.3 Void Galaxies 

In our SDSS DR7 galaxy catalog, there are 708,788 galaxies. 
In our M < —20.09 volume limited galaxy catalog, there are 
120,606 galaxies, with 8,046 of them falling inside voids, ap- 
proximately 7%. There are 79,947 (11%) void gala^xies that 
lie in void regions from the magnitude limited catalog with 
z < 0.107. Properties of these void galaxies will be discussed 
in a later paper. 
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Figure 5. The average radial density profile of all 1,054 voids in 
the void catalog as a function of the effective radius. 
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Figure 6. Radial density profi le (spherical annulus) as predicte d 
by linear gravitation theory llSheth fc van de Weveaer^ l2004l . 
The different curves correspond to different epochs of evolution, 
with the tallest peak representing 2 = 0. 



5 TESTS: VOLUME LIMITED CUTS 

In this section, we study the eflFect of changing the abso- 
lute magnitude cut on the voids found by VoidFinder. For 
absolute magnitudes brighter than Mr — —20.09, we use 
the same redshift cut while eliminating galaxies that fall 
under the absolute magnitude cut of —20.2, —20.3... — 20.6. 
For absolute magnitudes dimmer than Mr = —20, we use 
a redshift cut of a = 0.087, which corresponds to a limit- 
ing absolute magnitude of —19.5, and apply VoidFinder to 
samples with magnitude limits of —19.5, —19.6... — 20.1. It 
can be seen that as we slightly shift the absolute magnitude 
limits the void distribution remains similar, although there 
are trends that voids generally grow in size with brighter 
absolute magnitude cuts and voids get smaller in size with 
dimmer absolute magnitude cuts, as expected for changes in 
the sampling density of galaxies. We find qualitatively dif- 
ferent behavior as we examine extremely different samples 
(L*±0.5 magnitude), where we start to observe the effects 
of merging and splitting of voids. 

Figure [7] shows that the void regions found by 
VoidFinder are consistent for almost all large voids. The only 
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Figure 7. Overlap fraction for galaxy samples with redshift cut 
z = 0.107 (top), and z = 0.087 (bottom) with magnitude given in 
the figure compared to the void catalog sample (M;;^ = —20.09). 
The y-axis shows the fraction of the void volume that is also con- 
sidered void in the main sample as a function of the void volume. 
It can be seen that the large significant voids are consistently 
identified regardless of the volume limited cut. 



significant discrepancy arises from smaller voids that are in- 
troduced in sparser samples of the data. Figure [8] shows that 
the radial density profiles still show the "bucket shaped" fea- 
ture. 

Thus, the voids we find are not very sensitive to the 
absolute magnitude cut nor to the volume of our sample. 
SDSS DR7 provides a sufficiently contiguous three dimen- 
sional volume for void finding purposes. These voids found 
by VoidFinder should be considered significant large scale 
underdensities. 



6 TESTS: MOCK DATA 

We test the void finding algorithm on a set of mock galaxy 
catalogs to analyze the effects of the boundary conditions 
as imposed by SDSS, the effectiveness of studying large 
scale 3D structure in the finite volume of SDSS, as well as 
to test A-CDM predictions of the properti es of voids. The 
mock catalog used is a SPH halo model (|Skibba fc ShethI 
l2009h enclosed in a cube with sides 480 Mpc. The lu- 
minosity function and luminosity weighted correlation func- 
tions of the mock c atalogue are fit to SDSS as described by 
Skibba et al. Ill2006() , using halo occupation constraints from 



Zheng fc Weinberd l|2007l ) . 

We test VoidFinder on the mock sample using two dif- 
ferent methods. First, the SDSS mask is applied to the mock 
sample so that the geometries of the samples are the same; 
the results of this should mimic that of SDSS DR7. The 
SDSS geometry mock catalog contains 98,186 galaxies cov- 
ering a volume of 2.2 x 10^ Mpc"^ in the volume limited 
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Figure 8. Radial density profile (enclosed volume) for galaxy 
samples with -20.6 < Mr < -20.1 (top), and -20.1 < Mr < 
— 19.6 (bottom). The only difference between the profiles is the 
height of the peak at the edge of the voids. This is due to the 
different number density of galaxies in the sample used to deter- 
mine voids. The bucket shap ed behavior at the wal ls of the voids 
is consistent with ISheth &: van de Weveaer^ l|2004l 'l in Figure [6l 



Figure 9. Comparison of radial density profiles of voids in SDSS 
DR7 and simulations (top panel shows enclosed density, bot- 
tom panel shows density in spherical shells). The density profiles 
within the voids are nearly identical in all cases. The simulations 
show a slight tendency toward larger density just outside the void 
boundary. 



sample. Second, a cube is selected with volume similar but 
greater than the SDSS geometry sample. The cube geometry 
mock catalog contains 119,076 galaxies covering a volume of 
2.7 X lO^Mpc^ 
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6.1 Mock Results 

There are 1,006 voids and 6,228 void galaxies in the SDSS 
geometry mock catalog. There are 1,246 voids and 7,881 void 
galaxies in the cube geometry mock catalog. The void vol- 
ume fraction in the SDSS volume cut is 66.5%, and 69.3% in 
the cube volume cut. We observe that the geometry of SDSS 
plays a role in determining the overall void volume fraction, 
and if a SDSS geometry is considered in a mock sample, the 
volume fraction (66.5%) is approximately the same as the 
observed SDSS void volume fraction (62%). The effective ra- 
dius histogram of voids found in the mock samples in Figure 
1101 shows no significant changes in the sizes of voids found 
in the mock samples. The radial density profile in Figure [2] 
shows that the interiors of the voids are similarly empty as 
well. The void size and density profile results of the mock 
samples agree with observational data. However, there does 
seem to be a difference in the number of void galaxies found 
by VoidFinder which will be discussed in a separate paper. 



Figure 10. Distribution of void filling factor as a function of 
effective radius for voids found in SDSS DR7 and simulations. 
The distribution of void sizes found in the mock catalogs are 
nearly identical to those found in SDSS. The same size voids fill 
most of the volume. 
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7 VOID CATALOG RELEASE 

We have made this void catalog pubhcly available for future 
studies of voids. Included in the catalog are three separate 
interpretations of void regions. The first catalog consists of 
the maximal spheres of each unique void region. This is the 
largest hole in each void region in the shape of a sphere. This 
catalog is particularly useful for studying vast spherical un- 
derdense regions of the Universe. These spheres often depict 
the most underdense regions and galaxies near the centers 
of these voids are living in the most underdense large scale 
environments. The second catalog consists of all the possi- 
bly overlapping holes identified by VoidFinder. The merging 
of the holes forms each unique void region. This catalog is 
useful for identifying the entire void distribution of the Uni- 
verse. All of the volume enclosed by these holes lies in void 
regions and all galaxies contained are considered void galax- 
ies. The third catalog consists of the location and effective 
size of each unique void region. This catalog is useful for 
identifying overall void statistics in the Universe. Study of 
large scale structure as well as void volume distributions can 
be calculated from this catalog. Along with the three cat- 
alogs is the catalog of void galaxies. We have identified all 
galaxies with spectra that lie within the void regions iden- 
tified by VoidFinder. These catalogs can be downloaded for 
usfl 

We have now identified the largest and most compre- 
hensive void catalog from the largest spectroscopic data set 
available. Previous studies of voids from earlier data releases 
of SDSS, and other surveys including 2dF Galaxy Redshift 
Survey all lack the combination of completeness, depth, and 
contiguous sky provided in SDSS DR7. There is no longer an 
issue with survey boundaries restricting the volume of study 
for finding large voids. As there are currently no plans for 
a large spectroscopic survey of L* galaxies, this will be the 
most comprehensive data set for years to come. 



8 SUMMARY 

We studied the distribution of cosmic voids and void galaxies 
using Sloan Digital Sky Survey data release 7 using an abso- 
lute magnitude cut of Mr < —20.09. Using the VoidFinder 
algorithm as described bv lHovle fc"Vogele v ( 20o3), we iden- 
tify 1054 statistically significant voids in the northern galac- 
tic hemisphere greater than 10 Mpc in radius, covering 
62% of the volume. There are 8,046 galaxies brighter than 
Mr = —20.09 that lie within the voids, accounting for ap- 
proximately 6% of the galaxies, and 79,947 void galaxies 
(11.3%) with rrir < 17.6. The largest void is just over 30 
Mpc in effective radii. The median effective radius is 17 
Mpc. Voids of size r^ff ~ 20h~^ Mpc dominate the void 
volume. The voids are found to be significantly underdense, 
with 5 < —0.85 near the edges of the voids. We tested the 
sensitivity of the void finding algorithm to changes in the ab- 
solute magnitude cut within the range —19.6 > Mr > —20.6. 
The resulting void regions are largely similar with slight dif- 
ferences only near the edges of the void regions. The radial 
density profiles of the voids are found to be similar to predic- 
tions of dynamically distinct underdensities in gravitational 

1 www.physics.drexel.edu/~pan/ 



theory. We compared the results of VoidFinder on SDSS 
DR7 to mock catalogs generated from a SPH halo model 
simulation as well as other A-CDM simulations and found 
similar results, ruling out inconsistencies resulting from se- 
lection bias and survey geometry. 
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